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We study the discovery potential for detecting new neutral heavy Majorana leptons as suggested 
by some extentions of the Standard Model in recently proposed electron-proton colliders. Since 
1998-1999 the option of an electron-proton collider for the Very Large Hadron Collider at Fermilab 
operating with the proton booster has been considered. We consider the reaction e~p — > e + + jets 
and present estimates for the signal and Standard Model background including hadronization. 
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I. INTRODUCTION 

The recent experimental results on solar and atmo- 
spheric neutrino properties Q confirm that standard neu- 
trinos have small masses and oscillate. However, the na- 
ture of those massive neutrinos - if they are Dirac or 
Majorana particles - is still an open question. This is 
connected with lepton number conservation and to the 
general symmetries of some extention of the standard 
model. The smallness of neutrino masses is generally 
understood as a consequence of the see-saw mechanism. 
This brings the question of the possibility of new heavy 
neutrino states. So far none of these new states was ex- 
perimentally observed ]|, [| with masses up to Mjy ~ 100 
GeV. For higher masses there are many suggestions of ex- 
perimental possibilities in the next high energy hadron- 
hadron colliders [|[ in electron-positron linear accel- 
erators @, @, ||, in neutrinoless double-beta decay 
Another experimental possibility is found on electron- 
proton colliders with a positron in the final state [flo|| . 
This means a very clear experimental signature that vi- 
olates the electron leptonic number by two units. In this 
paper we investigate the discovery potential for recently 
proposed electron- proton colliders for detecting new neu- 
tral heavy Majorana leptons as suggested by extentions 
of the Standard Model. Since 1998-1999 the option of an 
electron-proton collider for the Very Large Hadron Col- 
lider at Fermilab operating with the proton booster has 
been considered. We study the reaction e~p — ► e + + 
jets and present estimates for the signal and Standard 
Model background. 



II. THE VERY HIGH ENERGY 
ELECTRON-PROTON COLLIDER 

There is a growing work on the options for a future 
Very Large Hadron Collider (VLHC)Oj. The low mag- 



netic field (approximately 2 T) collider design in a tunnel 
can also be used for a future high magnetic field (ap- 
proximately 10 T) collider. These designs could also be 
used for electron-proton colliders ]l2| [l|. In the low- field 
version electron-proton collisions would be produced by 
an 80 GeV electron ring and a 3 TeV proton booster 
(Va = 980 GeV), with a luminosity £= 2600 pb^/yr, 
while in the high-field option one could reach a center- 
of-mass energy ^/s = 6320 GeV and a luminosity L — 
1400 pb _1 /yr with the collision of a 200 GeV electron 
beam with a 50 TeV proton beam. These two options 
are considered in the present work. 



III. THE MODEL 

New heavy leptons are present in many possible exten- 
tions of the standard model. In the superstring inspired 
scenario 



14 we have new doublet vector leptons. In mir- 
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ror models | la] there are new fermions with right handed 
doublets and left handed singlets. Some see-saw exten- 
tions of the standard model include new weak-singlet 
fermions fUjfl . 

The general mixing of light and heavy Majorana neu- 
trinos has new mixing parameters as well as new pos- 
sible CP violating phases (lj]. It is well known that in 
neutrinoless double beta decay the theoretical estimates 
depend strongly on mixing angles and CP phases |l8[ . 
As a consequence, the comparison between the experi- 
mental limits and theory is model dependent. This is 
also the case for other processes. In this paper we are 
interested in the study of the discovery potential of new 
electron-proton colliders for new heavy Majorana neutri- 
nos. From the experimental point of view, instead of start 
comparisons of very particular models with a great num- 
ber of parameters and hypothesis one must start looking 
for the lightest of the new possible Majorana neutrinos, 
coupled to the electron family. This reduces the number 
of unknown parameters to two; the heavy neutrino mass 
and one new mixing angle, that should be determined 
exp er iment ally. 

In this paper we will adopt a general parametrisation 
that applies to all the above models. We will suppose 
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TABLE I: Branching ratios for N decays in the first family. 



Decays 


M N = 100 GeV 


M N = 1500 GeV 


N -> e+W" 


44% 


33% 


N -> e~W + 


44% 


33% 


N -f u R Z 


5,8% 


17% 


N -> i/ e Z 


5,8% 


17% 



that mixing angles for heavy-to-light neutrinos and new 
heavy neutrino masses are independent parameters p9| . 
The light neutrinos couplings to the neutral Z are given 
by gv,A — 9v A A — s * 7l2 ^mix/^' Taking the experimental 
values for gv,A and the Standard Model predictions we 
obtain a small upper bound for # m j x . A recent estimate 
§ gives sin 2 6> mix < 0.0052 with 95% C.L. This limit 
value is used throughout this paper for all curves and 
distributions. The decay of the Majorana neutrino into 
muons and taus are also possible. In ref. j|] it was found 
that for the muon sector the experimental mixing angle 
limit with 95% C.L. is much smaller than for the electron 
sector. 

Most extended models predict new fermions and new 
gauge vector bosons. Since we presently have no signal 
for new interactions, we will make the hypothesis that the 
new heavy neutrino states behave as SUl{2) ® Uy (1) ba- 
sic representations and interact mostly with the standard 
model gauge bosons. We can resume the new particle in- 
teractions in the neutral and charged current lagrangians: 

C nc = --—sinQmixZ^N^ (1 - 75) if} v + h.c. (1) 
and 

£cc = --^j^sinQmixW^NY (1 - 75) ike + h.c. (2) 

and similar terms for the other leptonic families. In the 
Majorana neutrino case, each completely neutral heavy 
lepton is coupled to all charged leptons but we are con- 
sidering only the lower mass state. 

The decay modes for these leptons, in the Majorana 
case j| must include both signatures N — ► £ T W ± and 
N — ► vt (i>g) Z, with £ = e, /i, t. In this work we have 
considered the Majorana neutrino decay only to the first 
family since for the other families we must add more mix- 
ing parameters. The Majorana neutrino branching ratios 
are shown in Table I and explicit formulas for the widths 
are given in ref. [po[ 



IV. RESULTS 

In Fig. 1 we show the structure of the Feynman dia- 
grams that contribute for the signal. We have included 



in our results all the Majorana off-shell and on-shell con- 
tributions. The calculations for the signal were done us- 
ing the high energy package CompHep plj with CTEQ95 
structure functions, after implementation of the extended 
model. We call attention to the fact that the N width is 
very narrow for lower heavy neutrino masses and increase 
for higher masses ( see Table II) . Thus propagator effects 
must be handled with care. 

The cross sections for on and off-shell heavy Majo- 
rana neutrino production as a function of its mass are 
displayed in Fig. 2 for y/s — 980 GeV and Fig. 3 for 
a/s = 6320 GeV. In all our results we have done first 
general detector cuts for the final positron E e + > 5 GeV 
and |77 e + 1 < 3.5, where rj e + is the pseudo rapidity defined 
relative to the initial proton. 

In order to reproduce the signal and background as 
close as possible to experimental data we have employed 
the package PYTHIA 6.1 [0 to generate the quark 
hadronization and decay from the Majorana neutrino 
production as well as the corresponding standard model 
background. The signal events were generated first at 
parton level with the CompHep package as above men- 
tioned and treated as external process to the PYTHIA 
program. In this way we have simulated the W, Z, quark 
and hadron decays. The Standard Model background 
was estimated leaving open all the possible channels in 
the PYTHIA program. The positrons come mainly from 
charm decay. In both cases we did not allow pions to 
decay. 

In the final visible particles we have imposed the fol- 
lowing criteria to select events: 

a) For \rj e + | < 3.5 and E e + > 5 GeV we must have the 
number of e + > 0, 

b) For \r)hi\ < 3.5 and E] li > 5 GeV we must have the 
number of hadrons > 5. 

With these two conditions, more than 99 % of the back- 
ground events were eliminated for both energies. 

In order to reconstruct the Majorana neutrino invari- 
ant mass we have employed the algorithm cluster PY- 
CELL from the PYTHIA package, with 60 azimuthal 
divisions and 26 divisions in the electron beam direc- 
tion. For y / s= 980 GeV we have applied the condition 
0.7 < V( A? ?) 2 + (^W whilc for 6320 GeV this 

condition was changed to 0.3 < \J (Ar/) 2 + (A0) 2 with 
1 77 1 < 3.5 for both energies. 

After these conditions, we have further demanded that: 

c) Number of clusters > 1, 

d) Number of clusters with 1 track = 1, 

e) Cluster invariant mass in the range (70, 90) GeV. 

This last requirement was done to identify the W com- 
ing from the heavy Majorana neutrino decay. After 
these additional cuts, the background estimated for 650 
k events for y/s ~ 1 TeV and 2,4 M events for y/s ~ 6 
TeV was reduced to practically zero. This point was ex- 
plicitly verified by generating 10 times more events than 
expected by the luminosity. For the signal, the same cuts 
imply the results shown in Table II. We also show in the 
last column of Table II. the significance S/y/S + B for 
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TABLE II: Number of events and significance for heavy Majorana neutrino production. 



yfs (GeV) 


M N (GeV) 


T N (GeV) 


Events after basic cuts 


Events after all cuts(%) 


S/y/S + B 


980 


100 


1.13 10 _;i 


516 


65 (12%) 


8 


300 


1.38 lO^ 1 


83 


32 (38%) 


5.6 


6320 


500 


6.50 10" 1 


334 


153 (45%) 


12.4 


1500 


17.6 10" 


80 


24 (30%) 
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each case, with " S" the number of signal events and " B" 
the number of background events. 

If we choose the clusters obtained with the surviving 
events then we can reconstruct the Majorana neutrino 
invariant mass as shown in Figs. 4 and 5. 

V. CONCLUSIONS 

The present work shows that heavy Majorana neutri- 
nos could be produced in the electron-proton option for 
the VLHC. The very clear signature of a positron in the 
final state allows to separate efficiently the signal events 



from the background and to reconstruct the Majorana 
neutrino invariant mass. In the yfs ~ 1 TeV option 
we can reach Majorana masses up to 300 GeV. For the 
yfs ~ 6 TeV option this value could be as high as 1500 
GeV. This is a conservative estimate since we are taking 
as limiting value S/y/S + B > 5. 
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Figure Captions 



1. General Feynman graphs for single heavy Majorana 



contribution to e p 



W - 



2. Cross section for single heavy Majorana neutrino 
production at yfs — 980 GeV using basic cuts 
E e + > 5 GeV and \r) e + \ < 3.5. 

3. Same as figure 2 for yfs — 6320 GeV. 

4. Reconstruction of the invariant Majorana mass dis- 
tribution for yfs = 980 GeV after hadronization 
using cuts and PYCELL algorithm. 

5. Same as figure 4 for yfs = 6320 GeV. 



